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ABSTRACT 

Repeated spectroscopic observations of stars in the Radial Velocity Experiment (RAVE) database 
are used to identify and examine single-lined binary (SBl) candidates. The RAVE latest internal 
database (VDR3) includes radial velocities, atmospheric and other parameters for approximately 
quarter million of different stars with little less than 300,000 observations. In the sample of ^20,000 
stars observed more than once, 1333 stars with variable radial velocities were identified. Most of 
them are believed to be SBl candidates. The fraction of SBl candidates among stars with several 
observations is between 10% and 15% which is the lower limit for binarity among RAVE stars. Due 
to the distribution of time spans between the re-observation that is biased towards relatively short 
timescales (days to weeks), the periods of the identified SBl candidates are most likely in the same 
range. Because of the RAVE's narrow magnitude range most of the dwarf candidates belong to the 
thin Galactic disk while the giants are part of the thick disk with distances extending to up to a 
few kpc. The comparison of the list of SBl candidates to the VSX catalog of variable stars yielded 
several pulsating variables among the giant population with the radial velocity variations of up to 
few tens of kms~^. There are 26 matches between the catalog of spectroscopic binary orbits (5^9 ) 
and the whole RAVE sample for which the given periastron time and the time of RAVE observation 
were close enough to yield a reliable comparison. RAVE measurements of radial velocities of known 
spectroscopic binaries are consistent with their published radial velocity curves. 
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1. INTRODUCTION 

Binary stars are not uncommon. Many studies search- 
ing for multiplicity among field and cluster stars of dif- 
ferent spectral types report that the fraction of binary 
stars in the observed sample is as high as 50% or more 
for certain spectral types. Multiplicity among dwarfs an d 
subdwarfs was st udied b y iDuquenno^ r fc MavoH dlflOlh: 



subawarls was stuaiea b y lUuciuennoy &: iviavoil lllblijlll: 
iFischer fc Marcvl ()1992[ ): iLadal (120061): iRaehayan et al.l 



(1201011 ■ ma s sive b inary stars bv iMason et al.l (|2009( ) 



Sana et all (120091 ) and binary stars in cluste rs by 
Abt fc WillmarthI ()1999[) : iSommariva et al.l (|2009[) . and 



others. Same findings are also supp orted by nu merical 
simulations of star cluster formation ( |Bate!l2QQ9( ). 

In the last few decades many new spectroscopic bi- 
naries were discovered, owing a large part to the suc- 
cessful CORAVEL and CfA speedometers and other 

|gal.matijevic@fmf.uni-lj.si| 



instruments. There were several surveys dedicated 
to the search of spectroscopic binaries (L atham et al] 
I2OOI iGrifRnI [200l iMermiUiod eFal] [20071 among oth- 
ers). The larger sca le Geneva- Copenhagen Survey 
(jNordstrom et al.l 120041 ) aiming at the F and G type 
dwarfs was also successful in identification of spectro- 
scopic binaries with a fraction of 19% in the observed 
sample. Many of the discovered binaries from various 
ca talogs along with their orbital parameters are compiled 
bv iPourbaix et aP (|2004[ ) in the 5*^9 catalog of spectro- 
scopic binary orbits. 

With the development of optical-fiber spectrographs 
that enable the observation of up to several hundred stars 
simultaneously, efficient large-scale spectroscopic surveys 
became possible. Among the largest such sky surveys 
with the focus on stellar objects are the Sloan Digital Sky 
Survey (SDSS) covering the nort hern sky and the Ra,- 
dial Velocity Experiment (RAVE; iSteinmetz et al.ll2006f l 
covering the southern sky with a higher spectral resolu- 
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tion but narrower wavelength range than the SDSS. Bi- 
nary star ca ndidates observed as a part of the SDSS are 
discussed in iPourbaix et all ()2005l ) . Double-hned spec- 
troscopic binary (SB2) candidates discovered in the sec- 
ond data release of the RAVE survey where published in 
iMatiievic et al.l ()2010t ) - Paper I hereafter. The analy- 
sis of single- lined spectroscopic binary (SBl) candidates 
identified in the RAVE's latest internal database (VDR3) 
is the subject of this paper. Section [2] gives an overview 
of radial velocity (RV) acquisition in the RAVE survey. 
Section [3] discusses the sample of RAVE stars that were 
observed multiple times and the procedure that was used 
to identify potential SBl objects. In Section^ a catalog 
of SBl candidates is presented. The concluding section 
summarizes the main results and discusses future work 
to use the discovered binary candidates for population 
analysis and to assess their influence on the Galactic RV 
distribution. 

2. RAVE RADIAL VELOCITIES 

RAVE is an ongoing multi-fiber spectroscopic survey 
based on observations with the UK Schmidt Telescope 
at the Australian Astronomical Observatory (AAO) with 
the goal of observing up to 1 million stars in the mag- 
nitude range between 9 < / < 12. The wavelength 
range of the spectra covers the near-infrared region 
AA8410 - 8795 with a resolving power of i? - 7500, typ- 
ically with a moderately high signal-to- noise ratio (S/N; 
mean value of 45 pixel"^). The RAVE spectral analysis 
pipeline is designed to derive accurate stellar RVs as well 
as atmospheric pa rameters and chemical composition 
(jBoeche et al."2010 D. So far^ two d a ta relea ses have been 
published (Stcinm etz et all l2006t iZwitte r ct al. 2008) 
and the third one was submitted (jSiebertet al.ll2011i) . 

The measurement of stellar RVs is one of the major 
goals of the RAVE survey. Along with the star's position 
on the sky, its proper motion and distance, an RV gives 
us a missing sixth kinematic component. Having full six- 
dimensional data for up to a million stars allows one to 
carry out detailed studies of structure and kinematics of 
the Galaxy. 

All data provided by the RAVE survey are acquired 
at the AAO 1.2m UK Schmidt telescope with the 6dF 
optical fiber system feeding the dedicated spectrograph. 
About a hundred optical fibers are positioned in the fo- 
cal plane of the telescope to route light from the same 
number of stars simultaneously to the slit vane mounted 
in the spectrograph enclosure. Additional fibers are used 
to observe the background sky. The slit vane stacks the 
endings of all fibers in a column, one over the other. Due 
to the close spacing of the fibers there is a chance of 
light cross-talk between them. Light from one fiber may 
be spilling into the adjacent one. This effect is minimized 
by observing stars in a narrow magnitude range and by 
data reduction pipeline which subtracts the cross-talk by 
an iterative procedure. Second order light which contam- 
inates spectra from the first RAVE data release is now 
blocked by a blue blocking filter placed in front of the col- 
limator. The light passing the collimator is dispersed us- 
ing a volume phase holographic transmission grating and 
finally recorded with a camera using a GGD detector. Be- 
cause the throughput of the different fibers may vary, the 
final S /N of two equally bright stars from the same field 
plate may be different. Typically, fields are observed in 
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FlG. 1. — Cumulative distribution of measured radial velocity 
errors. The dashed line represents the data from the VDR3 release 
with S/N > 20. The thinner black line shows a subset that includes 
only stars that were confidently classified as cold (T^g < 7500 K) 
normal single stars. The dotted line represents only the already 
published data from the second data release. Thicker light gray line 
shows the distribution of applied zero-point corrections (to correct 
the temperature variations, see the text) and the thicker dark gray 
line shows the average discrepancy between RVs and their weighted 
average in terms of RV errors (Equation jlj) for stars with multiple 
observations. Data from the first release plagued by second order 
light contamination are excluded in all shown sets. Dash-dotted 
lines mark the 50th, 68th and 95th percentiles. 

the following manner. Before and after five consecutive 
field exposures lasting 10 minutes each, two calibrating 
lamp exposures are usually made. Slight temperature 
variations in the spectrograph enclosure may shift the 
wavelength calibration so each plate has a few dedicated 
sky fibers that (together with sky velocities derived from 
stellar spectra) serve as additional wavelength constraint 
and are also used for the subtra ction of the tellu ric lines 
present in stellar spectra (Stein metz et al.ll2006( ). 

The measurement of RVs is performed in a following 
way. First, a spectrum is continuum normalized us- 
ing a pair of cubic splines. Then, a four step iterative 
method is used to calculate its RV: (1) using a stan- 
dard cross-correlation procedure (jTonrv fc Davis! [l979( ) 
a first RV estimate is calculated. This step requires a 
library of template synthetic spectra against which a 
correlation is calcul ated. A subset of 40 spectra from 
iMunari et ahl (|2005[) covering the whole {Teg, log pa- 
rameter space is used to get a rough estimate with a 
precision of ~ 5 km s~^ for spectra with strong Ca II 
lines. (2) The spectra are shifted relative to the rest 
frame according to the RV estimate just calculated. (3) 
A best-matching template is constructed fr o m a f ull li- 
brary of template spectra (see IZwitter et al.l ()2008D for a 
more detailed description of this procedure). (4) With 
a best-matching template a new RV is calculated, and 
finally this measurement is corr ected for a possible zero- 
point shift (jSiebert et al.ll20Tll ) that is produced by the 
temperature variations in the spectrograph components. 
Errors of the RV measurements are estimated by fitting a 
parabolic curve to the peak of the cross-correlation func- 
tion. 

The VDR3 database consists of 344,924 entries includ- 
ing the already published data from the first and second 
releases. Out of those there are 295,618 database en- 
tries with S/N > 20 and further exclusion of the second- 
order-light-plagued data from the first release leaves 
279,120 observations of 249,980 different stars. This sub- 
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set serves as the basis of this paper. Repeated observa- 
tions were identified by matching Two Micron All Sky 
Survey (2MASS) identifiers. For a small set of stars 
with missing 2MASS identifiers we compared their co- 
ordinates and matched stars that were not more than 5" 
apart. All matching cases were visually checked on the 
DSS plates to ensure the catalog entries indeed belong to 
the same star. We decided to exclude the data from the 
first release because the classification of those spectra is 
unreliable due to strong continuum variations and their 
RV measurements are usually significantly less accurate 
than the RVs of spectra recorded later. Cumulative dis- 
tributions of RV errors for different data sets are shown 
in Figure [TJ The median error of the selected subset 
stands at just over 1km s~^ and RVs of 95% of all cold 
single stars are measured to better than 2 km s~^. There 
is a notable improvement in RV accuracy compared to 
the second data release. The same figure also shows the 
distribution of applied zero-point corrections. In roughly 
two thirds of all cases this correction is smaller than the 
error of the measured RV, although in some cases the 
correction can be an order of magnitude larger than the 
RV errors. 
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3. REPEATED OBSERVATIONS 

In order to identify any stars with variable RVs, more 
than one observation of the same stars is needed. The 
RAVE survey is primarily focused on the acquisition of as 
many different targets as possible. Re-observations of the 
same stars are done for stability and quality check pur- 
poses and also to allow for the determination of the frac- 
tion of stars with variable RVs. Nevertheless, at about 
9% of all observations the number of stars that were ob- 
served more than once is still high (Figure[2|). Altogether, 
21,730 stars were observed more than once in the selected 
subsample. Most of the repeats were observed in the fol- 
lowing few days with a strong peak at a day after the 
previous observation. 

The stability of RVs is shown in Figure [1] The thick 
dark gray line shows the distribution of average differ- 
ences between the RV measurements and their weighted 
average (RV) in terms of individual errors, 
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for all repeatedly observed normal single stars with 
S/N > 20 (same set as used for the thin solid line in 
Figure [T]). N indicates the number of observations for 
a given star. A long tail of the distribution toward the 
larger values of A is indeed not an observational error but 
indicates the presence of the stars with variable RVs. 

To establish a quantitative criterion for RV variabil- 
ity we defined the following function. We denoted two 
RVs and their errors measured for the same star at two 
different times as (RVi,cri) and (RV2,(T2)- The squares 
of the RV errors u1 can be treated as variances of the 
Gaussian distributions with the mean values RV^. If we 
would randomly pick two samples from each of these dis- 
tributions, the probability that the pick from the second 
one is greater than the pick from the first one is equal to 



Fig. 2. — Upper diagram shows the number of objects observed 
n times along with the fraction that this number represents in 
the whole sample of 279,120 stars. The lower diagram shows the 
distribution of the time span between the consecutive observations 
of the same objects. 
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The double integral can be simplified by introducing a 
new variable u = [x — RVi)/-\/2cri and evaluating the 
integration over x, 



P{2 > 1) 



1 



1 + erf 



y-RVi 
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(3) 



(B-RVa) 

e dy, 



where erf is the standard error function. By using the 
following identity. 



eTf{cx+d)dx 



/TT I ad — he . 
-erf ( ^== 1 , (4) 



we can write a simple expression for the probability in 
question, 



^(2 > 1) = ^ 




(5) 



If the difference in the numerator of the error function is 
zero (both RVs are the same), the probability is equal to 
1/2 and goes toward 1 for pairs with very different RVs 
and comparably small as. If RVi is greater than RV2, 
than the probability P will be less than 1/2. Since we 
are only interested if the RV changed between measure- 
ments, we can always label the greater of both RVs with 
index 2 without the loss of generality. The nature of er- 
ror function makes it hard to tell the difference between 
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Fig. 3. — piog from Eq. |6]as a function of difference RV2 — RVi 
in units of ai. Different line styles correspond to different values 
of (T2 , also in units of ci . 

two pairs with relatively large RV differences because the 
value of the error function will always be very close to 
1 in such cases. It is more convenient to define a new 
function involving the logarithm of the function P, 



piog = -logio(l - P). 



(6) 



The calculation of the logarithm in this function where 
the difference 1 — P 10~^^ is troublesome due to dou- 
ble precision floating point limitations. To prevent any 
numerical inconsistencies the upper limit of the piog cri- 
terion was set to 14 in the subsequent calculations. The 
function piog is shown in Fig. [3] for several different error 
ratios. For example, if both errors are the same (ui — 02) 
and the RVs are roughly 4cti apart, the value of the func- 
tion piog will be a little less than three, corresponding to 
the value of P = 0.998 so the large majority of the values 
from the second distribution are greater than the values 
from the first. The probability that the radial veloc- 
ity changed from one observation to another is relatively 
high. On the other hand, if piog < 2, the variability is 
questionable and if piog < 1, the radial velocity variabil- 
ity is insignificant. In compar ison, the lower l i mit of the 
variability criterion given in IPourbaix etahl (pOOl for 
stars observed twice and with equal RV errors for both 
measurements is piog — 2.87. 

4. CATALOG OF SBl CANDIDATES 

We have examined a sample of 21,730 objects that were 
observed at least twice. We excluded all objects whose 
spectra are morphologically different from spectra of sin- 
gle stars since the described RV extraction method does 
not give reliable results in such cases, including double- 
lined spectroscopic binaries. All such objects were iden- 
tified using the classification method described in Pa- 
per I. The method systematically examines the proper- 
ties of the cross-correlation function between each ob- 
served spectrum and a pre-calculated synthetic template 
and according to those properties groups similar spec- 
tra together. Among the excluded classes were spec- 
tra with observational or reduction errors (problematic 
continuum normalization, for example), hot stars whose 
Paschen series lines are too wide for precise RV measure- 
ment and intrinsically peculiar stars (stars with active 
chromospheres, emission line stars, etc.). The selection 
reduced the number of stars in the sample to 20,027. 

We calculated the piog variability criterion for all stars 
from this sample using the RV measurements and their 
errors as given in the RAVE database. In cases where 
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FlG. 4. — Bottom diagram shows the ratio between the number 
of SBl candidates vs. the number of the rest of repeatedly ob- 
served stars as a function of zero-point shift. All histograms are 
normalized so that the first bin equals 1. Gray lines show the initial 
calculation w here the zero-point uncertainty is not included (fe = 
in Equations I7I8I I. whereas black lines show the same distribution 
for k = 0.5. SBl candidates in both cases are considered all pairs 
with piog > 3 (full lines) and pi^g > 4 (dashed-lines). The upper 
diagram shows the distribution of zero-point shifts for the pairs 
of observations with the highest variability criterion for all stars 
observed more than once. 

more than two observations were available, the pair with 
the highest value of the criterion was considered. The 
inspection of a few stars with the largest number of re- 
peated observations revealed that in some cases only one 
or two of the RV measurements were not consistent with 
the averaged RV. We note that some of the observations 
responsible for the RV variability have large values of the 
zero-point correction which is usually very small (Section 
[21). The comparison of the number of SBl candidates 
with zero-point shift in some interval to the rest of the 
repeatedly observed stars in the same interval showed 
that there was a significant overdensity of SBl candi- 
dates in the region where at least one of the two RV 
measurements in the pair was corrected for a relatively 
high amount (Figure 2]). Such spurious detections would 
contaminate the statistics of discovered SBl candidates. 
To account for this systematic error we recalculated the 
variability criterion with the adjusted error estimate: 



RV 



'ZP' 



(7) 



where 



CTzp ^ k* RVzp, 



(8) 

so the error in zero-point shift is proportional to the shift 
in this approximation. Equation[7]is only valid if RV and 
zero-point error estimates are uncorrelated. This was 
checked and indeed holds true. The coefficient k = 0.5 
was selected so that it fully eliminates the SBl overden- 
sity at large zero-point shifts. It was also intentionally 
set high enough so that the ratio of SBl candidates con- 
firmed with the adjusted error estimate compared to the 
ones before this test is a decreasing function of the value 
of the zero-point correction. One expects that because 
high values of zero-point shifts are the least trustworthy 
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TABLE 1 

Number of SBl Candidates for. Different Values 

OF piog 
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and include the highest fraction of potential false posi- 
tives. The overall number of SBl candidates is not much 
affected by the selected value of k, since the majority of 
all RVs have relatively small zero-point corrections (in 
only 5% of the cases they exceed 3 km s~^ and anyway 
are smaller than the RV errors, see Figure [T]). 

Altogether, 1333 or 6.6% of the stars have piog > 3 
and were identified as potential SBl candidates. The 
summary of the number of discovered candidates versus 
the number of observations is given in Table [TJ Since the 
limiting value of piog = 3 is arbitrary, p\og = 2,4 are given 
for comparison. It is evident that for piog = 2 the number 
of SBl candidates becomes unrealistically high, making 
the piog = 3 a plausible lower limit for the variability cri- 
terion. The efficiency of the detection is ~ 6% for twice 
observed star but it grows steadily toward ~ 15% for 
stars observed five or six times. There are too few can- 
didates with more observations to draw any conclusions. 
The number of identified candidates gives only a lower 
limit of spectroscopic binaries in the observed sample and 
selection criteria were indeed set restrictively to avoid as 
much potential false positives as possible and also to ac- 
count for the fact that the RV error estimates might not 
always be Gaussian. With the maximal time span be- 
tween re-observations at around 2000 days it cannot be 
expected that systems with significantly longer periods 
are detectable. More so, in long period systems the RV 
amplitudes become too small to be detectable and most 
of the SBl candidates were observed for the second time 
only days after the first observation. The distribution of 
RV differences between measurements is shown in Fig- 
ure [S] Most frequently the variations were small, around 
the value of 6 km s~^. The distribution falls quickly af- 
ter 10 km s~^ but extends all the way to over 160 km s^^. 
All single-lined objects with high RV shifts are particu- 
larly interesting for further investigation due to the pos- 
sibility of the presence of massive and faint objects. The 
distribution of RV shifts was also compared to the dis- 
tribution of separations between the components of the 
identified SB2 candidates (Figure [S]) . The data for ve- 
locity separations were taken from the list of preliminary 
solutions for 1040 SB2 objects found in the RAVE catalog 
(partially described in Paper I). The distribution shown 
on the plot was scaled compared to the SBl distribution 
so that the effective number of stars from which the selec- 
tion was done is the same for both types. The efficiency 
of SB2 detection becomes high at the point where the 
number of identified SBl objects gets low. 
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Fig. 5. — Distribution of maximum changes of radial velocity 
between measurements for the identified SBl candidates. The 
smaller diagram shows the same distribution in the logarithmic 
scale along with the preliminary normalized distribution of sepa- 
rations between components of identified SB2 candidates (dashed 
line). 

The {Toff, log 5} diagram of all SBl candidates and 
various distributions are shown in Figure [Bl The dis- 
tribution of effective temperature of stars observed mul- 
tiple times has two distinct peaks at around 4500 K for 
the red clump stars and 6000 K for the main sequence 
dwarfs, same as in the overall RAVE sample. The S/N 
of the re-observed stars is somewhat higher than in the 
general population. The reason is that slightly brighter 
stars were observed more frequently than the faint ones. 
There are more SBl dwarf candidates than giant candi- 
dates, which is also apparent from the distribution of ef- 
fective temperatures for the SBl candidate sample. The 
ratio of giants to dwarfs in the whole RAVE sample is 
around 57:43, while the same ratio for the SBl candi- 
dates is close to 50:50. This is an expected result. SBl 
candidates roughly fall into two groups: main sequence 
stars with masses ~ 1 — 1.2 Mq and red clump giants 
with masses slig htly larger than the first group (based 
on isochrones bv iMarigo et al.l ()2008[ )V The important 
difference between the two groups is approximately an 
order of magnitude larger radii of giant stars. This im- 
plies that the smallest binary orbits than can host giant 
stars must be larger than the orbits of main sequence 
stars due to Ro che lobe radius limits (see for example 
lEggletonl (|1983f )l. Consequently, the periods of SBl bi- 
naries that host giant stars will be longer and average 
RV shifts will be smaller and therefore harder to detect. 

The distribution of metallicity is relatively wide since 
RAVE observes field stars representing the general pop- 
ulation of thin and thick disks. The metallicity of SBl 
candidates does not differ from the rest of the stars mean- 
ing that SBls are also scattered throughout both disks. 
Due to rave's narrow magnitude range, dwarf candi- 
dates are mostly limited to the thin disk, while giants 
also reach into the Galaxy's thick disk (lower left panel 
of Figure [SI see also iZwitter et al.l (|2010[ ) for further dis- 
cussion on distances). 

Visual inspection of spectra revealed that 43 SBl can- 
didates could in fact be double-lined spectroscopic bina- 
ries if observed at higher resolutions or at more favorable 
phases (red markers in the upper right panel of Figure 15]) ■ 
The assumption is based on their slightly asymmetric and 
widened lines but the effect is too small for more confi- 
dent confirmation. Note that atmospheric parameters 
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Fig. 6. — Upper left panel: {Tefijlogg} diagram for a sample of identified SBl candidates. The variability criterion (piog) is represented 
with different shades of gray where darker tones and bigger markers correspond to greater chance of RV variability. The dashed line roughly 
divides giants mostly from the red clump (RC) region and main-sequence (MS) stars. The error bars represent typical uncertainties of 
both parameters. Upper right panel: photometrically variable stars from the VSX catalog (see the text) and suspecte d SB2 objects among 
the identified SBl candidates. Lower left panel: SBl candidates with color-coded distances from lZwitter et al.l 1)201 0) and tone values as 
in the previous diagram. Here, only the size of markers corresponds to the value of piog to avoid confusion. Isochroncs for [M/H] = —0.2 
spannmg 9.0-10.0 in log age with 0.2 steps bylMarig o et al., (,200S) and lGirardi et al.l 1I2OIOI '). Lower right panels: normalized distributions 
of effective temperature, metallicity, 2MASS Mj magnitude and S/N for non-peculiar single stars of the whole RAVE sample (blue), stars 
observed more than once (green) and SBl candidates (red). 



for such stars are generally incorrect since the pipeline 
always assumes that the observed spectrum belongs to a 
single star and models the binary spectrum with this as- 
sumption. Nevertheless, the general region (dwarf/giant) 
of the solution is correct since the most affected parame- 
ter is rotational velocity while the effective temperature 
and surface gravity are not very sensitive to the appar- 
ent slight line broadening that is produced by the shift. 
Most of these stars are close to the main sequence. 

We do not wish to attempt to infer binary orbits for 
those stars with the highest number of observations. 
While it is possible to derive orbital parameters for stars 
with six or more observations in theory, the precision 
of the calculated solution also depends on how well the 



measured points are distributed along the orbital phase, 
on the uncertainty of RVs and so on. We tried to fit the 
orbits for some stars, but none of the solutions turned 
out to be reliable. 

The catalog of all stars with variable RVs will be made 
publicly available through the CDS service. 

4.1. Comparison to other catalogs 

We compared the list of identified SBl candidates 
to the li st of known ph otometric variables in the VSX 
catalog (jWatsonl 120(361 Version 02-Jan-2011) and to 
the ninth catalog of spectroscopic binary orbits {Sb^', 
iPourbaix eraLll2004D . Within the VSX there are 36 stars 
matching our list. Out of those, 10 are pulsating vari- 
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Fig. 7. — Spectroscopic orbits from Sgg catalog and RAVE RVs (black mark ers) for each star. The p hase uncertainty is calculated from a 
given period error. The radial velocity data for HIP 5438 from lAndersen et al.' ( 1988), HIP 44124 from Goldberg et al? (1002), BD -21 3873 
from Smith et alj l(T997l ). BD -01 3022, H IP 43527 and CPD - 64 4333 from Udry et al. ( 1998), HIP 75718 and HIP 86400 from Tokovini^ 
||T991I ). HIP 59750 and HIP 55022 from ICarnev et al.l l[200lD . BD -01 3141 from Griffin (1994), HIP 65982 from Latham et al. (2002), 
CI* NGC 2477 HART 206 4 and CI* NGC 2477 H ART 80 1 7 from [Eigenbrod et al. (2004), HIP 42368 from Griffin (2004), HIP 64219 
from IDuouennov fc Mavoa ( IT9911) . HIP 50796 from ITorrej 1200^) . BD -18 5091 from iCarguillat Prieurl ( i2007i ). BD -20 3 310, BD - 
15 3197, CD -23 9811, CD -27 9 246, CD -27 9317 , CD -26 9873, CD -25 10022, BD -20 3997 and HIP 22761 from IGriffiiJ i200B) and 
CI* Ruprecht 46 MMU 2191 from lMermilliod et al.l ll2007l ). 



8 



Matijevic et al. 



ables, 14 eclipsing binaries and 20 unclassified variables. 
Their positions on the {Tgff, log diagram are shown in 
Figure [6] RV variations measured for pulsating variable s 
are generally lower than 30km s""'^ (lAerts et alJl2010l ). 
but can be as high as ^ 100 km s~^ for some types 
(|Gorvnva et aUfTOOSt ISanford I ri9491 In our case, the 
two leftmost green dots in the upper right panel of Fig- 
ure [6] are an RR Lyrae star (hotter) and a Cepheid with 
the RV amplitude of at least 21km s~^ and 46 km s~^, 
respectively. The rest of the pulsating variables mostly 
lie above the dwarf-giant border. Of course, there is a 
possibility that some of the rest of the identified SBl can- 
didates are pulsating variables rather than binary stars. 
Conversely, eclipsing binaries and potential SB2 candi- 
dates are generally (taking atmospheric parameter uncer- 
tainties into account) among main sequence dwarfs. This 
supports the assumption from Paper I an d also agrees 
with t he sample of eclipsing binaries from iTorres et al.l 
(|2010f) where the majority of all analyzed systems are 
pairs of main-sequence stars. 

The comparison of the whole RAVE sample to the 5*^9 
catalog yielded 56 matches but only three of those stars 
were observed more than once, two of them twice, and 
one four times. For the sake of testing the accuracy of 
RVs, single observations were compared to known orbits 
as well. We selected 26 orbital solutions for which the 
difference between the published periastron time and the 
date of the RAVE observation were closer than 100 peri- 
ods apart in order to avoid errors caused by period propa- 
gation. There are two exceptions with a difference of 343 
and 448 periods where the correspondence between the 
RAVE measurements and the RV curve was still accept- 
able. In other cases several thousand periods have passed 
between the given periastron time and the time of the 
RAVE measurement. Orbital solutions with RAVE mea- 
surements are shown in Figure [71 There are four known 
SB2 systems among the selected sample. RAVE spectra 
of HIP 5438 and HIP 44124 were recorded close to the 
half phase and so both RVs are near the systemic 7 veloc- 
ity. The RV amplitudes of both systems are large enough 
that both systems would be detected as SB2 if observed 
at more favorable phases. System BD -20 3310 was ob- 
served close to the apastron but again the measured RV 
is near the center-of-mass RV because spectra of both 
components contribute to the composite so the peak of 
the correlation function lies in the middle of both com- 
ponents. System CD -27 9246 was observed four times, 
but coincident ally at only two different phases. All four 
RVs correspond to one component of the system. Other 
RAVE RVs are close to the predicted values and in most 
cases within the error bars, giving additional reassurance 
to the error estimates. In three cases, BD -21 3873, BD - 
01 3141, and HIP 22761, the measured RVs lie very close 
to the 7 velocities of the system, very similar to the three 
cases of known SB2 systems. The reason for this seems 
to be different in each case. A poor match in the case 
of HIP 22761 (HD 31341) is most likely caused by con- 
fusion with another star due to a n erroneous e ntry in 
the catalog. The S'^g catalog cites iGrifhnl ()2006[ ) as the 
source, but the paper actually reports about the obser- 
vation of the star HD 31341B, also known as HIP 22766. 
The mismatch in BD -21 3873 could be caused by an 
improper phase determination of the two RAVE mea- 



surements, since the phase inaccuracy in this case is the 
largest of all examples. In the case of BD -01 3141, a pos- 
sible secondary component might only be observable at 
longer wavelengths because of its pote ntially low surface 
temperature. The reference paper (Griffinl Tl994[ ) states 
that all observations were performed on instruments ca- 
pable of reaching wavelengths up to '-^ 5200 A which is 
relatively far from RAVE's spectral domain so there is 
a chance that the additional component was not present 
there. 

5. DISCUSSION AND CONCLUSIONS 

We analyzed 20,027 stars observed multiple times by 
the RAVE spectroscopic survey in the ^ 6 year span be- 
tween 2004 and 2010 with RV and atmospheric parame- 
ter estimates provided by the RAVE parameter estima- 
tion pipeline. Prior to the RV analysis we classified all 
spectra and filtered out those with peculiar features (e.g., 
emission line objects, double-lined binaries, spectra with 
observational errors, etc.) whose measured parameters 
are not reliable, obtaining a list of only normal single 
stars. In order to detect stars with variable RVs, we 
defined a variability criterion (Equation ([5])), assuming 
that RV errors are Gaussian and taking the RV error and 
zero-point shift of individual measurements into account. 
In the observed sample, we identified 1333 objects with 
large enough changes in radial velocities to be identified 
as SBl candidates. As summarized in Table [1] the frac- 
tion of discovered SBl candidates is a function of the 
number of observations. Only 6 % of twice observed 
stars were detected as SBl candidates. This fraction 
grows with the number of observations and saturates at 
around 10%-15% for stars with five or six observations, 
depending on the selected lower limit of the variability 
criterion. This number represents the lower limit for the 
overall fraction of stars with variable RVs present in the 
RAVE sample. 

The distribution of maximal differences between the 
RVs for repeatedly observed stars (Figured]) has a strong 
peak just below 10 km s~^. At smaller velocities the num- 
ber of detected SBl candidates quickly vanishes. The 
reason is a restriction on the piog criterion that is re- 
quired to be greater than 3 for sufficiently confident 
variability confirmation. Even in cases with the low- 
est values of RV errors (Figure [1]) this criterion is not 
met if the two measurements are closer together than 
~ 3.5kms~^ (Figure 131) which explains the sharp falloff 
below 4kms~"'^. The same distribution was plotted for 
stars where piog > 4 and the falloff was even stronger 
there which gives additional confirmation to the source 
of this feature. It should be noted that zero-point shifts 
(Figure [IJ are in most cases well below the velocities 
where the detection becomes significant and are not the 
source of inaccuracies. On the other side of the distri- 
bution, there is an evident exponential falloff of detected 
candidates. This feature is expected to be real and not 
a selection effect. Systems with larger RV amplitudes 
most likely have shorter periods and therefore the de- 
tection efficiency for these cases is proportionally higher 
than at the lower velocity end. The distribution also ex- 
tends into the SB2 region for which it was shown (Paper 
I) that the detection efficiency is very high at separa- 
tions above 50kms"^. It is hard to infer the periods 
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of the SBl candidates since mostly only two observa- 
tions are available, but from the distribution of the time 
spans between observations it can be concluded that pe- 
riods are most likely not significantly longer and so fall 
into the short end o f the p eriod distri bution g i ven i n 
iDuquennov k Mavoi] (|1991l ) or iRagha van et alJ (|2010l ). 
A detailed population study of both types of binaries 
(SBl and SB2) and their influence on the general RV 
distribution of the RAVE stars will be a subject of the 
third paper regarding binaries in the RAVE survey. 

The {Toff, log g} diagram of SBl candidates derived 
from the spectroscopic parameters of stars gives a similar 
picture as RAVE's general population. The distribution 
of metallicities of SBl candidates is the same as the dis- 
tribution of metallicities of the general population which 
means that binaries are well mixed among the field stars. 
Similarly to the overall RAVE sample, the distribution of 
effective temperatures has two distinct peaks that corre- 
spond to main sequence stars and to the red clump giants 
with masses only slightly larger than the former group. 
There is a deficiency of SBl candidates from the giant 
group compared to the dwarfs. It can be explained by 
comparing the typical radii of stars in both groups. Radii 
of giants in the observed sample are roughly an order of 
magnitude larger than the radii of main sequence stars. 
Therefore, binary orbits in which giants can exist (due 
to Roche lobe limits) must be on average larger than 
the orbits of main sequence stars so the orbital periods 
are longer and consequently RV shift smaller and more 
unlikely to be detected. 

Comparison of the SBl candidate list with the VSX 
catalog of photometrically variable stars yielded several 
hits with some stars being pulsating variables rather than 
spectroscopic binaries. The RV amplitudes of some types 
of pulsating variables can be as high as few lOkms"^. 
There is some probability that some of the stars from our 
SBl candidate list are pulsating but it should be negligi- 
bly low. For example, the fraction of Cepheids and RR 
Lyrae stars (variables with high RV variations) in the 
Hipparcos catalog (Ferryman et al. 1997) is less than 0.4 
%. The selection of stars in the RAVE's input catalog 
is unbiased (with the exception of stellar magnitudes) so 
most of the observed stars are not in any unstable tran- 
siting phases of their lives. Among the photometrically 
variable stars are also several eclipsing binaries. Most 
of them lie close to the main sequence on the HR dia- 
gram. We also identified some potential SB2 candidates 
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in our SBl candidate list. They cannot be confirmed 
since they do not have a clear double-lined signature in 
their spectra, but time dependent asymmetries in spec- 
tral lines indicate that higher resolution spectra might 
be able to resolve two components. All of these objects 
also lie close to the main sequence which supports the 
argument that SB2 objects mostly consist of two main 
sequence stars and pairs of equally bright giants forming 
an SB2 system are rare. 

There are 26 stars in the catalog of spectroscopic bi- 
nary orbits (5^9) for which RAVE has at least one ob- 
servation and the difference between the given time of 
periastron and the RAVE observation time is shorter 
than 100 epochs to ensure that the propagated period 
error is small enough. In most cases the agreement be- 
tween the calculated RV curves and RAVE RVs is within 
the error estimates which confirms the precision of the 
RVs. Among the selection are four known SB2 objects. 
In three of those cases, the RAVE RV is close to the 
7 velocity of the system which is expected because the 
peak of the correlation function that determines the RV is 
between both Doppler-shifted components. Conversely, 
in three other cases that are cataloged as SBl objects, 
RAVE radial velocities are close to the systemic velocity 
and not near the calculated RV curve. This might impli- 
cate the presence of a still undiscovered bright secondary 
component. 
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